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Abstract 

We consider a noncommutative description of graphene. This description consists of a Dirac 
equation for massless Dirac fermions plus noncommutative corrections, which are treated in the 
presence of an external magnetic field. We argue that, being a two-dimensional Dirac system, 
graphene is particularly interesting to test noncommutativity. We find that momentum noncom- 
mutativity affects the energy levels of graphene, but that it does not entail any kind of correction 
to the Hall conductivity. 
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I. INTRODUCTION 



Graphene is a two-dimensional configuration of carbon atoms organized in a liexagonal 
lioneycomb structure [l-^. Often, a crystal lattice is a Bavrais lattice, that is, an infinite 
array of discrete points with an organization and orientation that appears exactly the same, 
from whatever point the array is viewed. However, the graphene hexagonal lattice is non- 
Bravais, as only the next-to-nearest neighbor points appear with the same organization 
and orientation. In the case of graphene one has two triangular Bravais lattices, A and 
B which together form the non-Bravais graphene lattice, and the difference between them 
is a rotation of vr. The hexagonal lattice belongs to the class of bipartite lattices, and 
so one can say that graphene is a bipartite non-Bravais lattice with two carbon atoms 
per unit cell. For the two sub-lattices one has the same primitive vectors, which depend 
explicitly on the distance between the two lattice points. With these primitive vectors one 
can characterize any space point as a linear combination of them. In the momentum space, 
one can obtain the reciprocal primitive vectors and the corresponding Brillouin zones, i.e. 
a uniquely defined primitive cell in the reciprocal space. The first Brillouin zone forms a 
hexagon, which is rotated by 7r/12 compared to the hexagonal structure in position space 
and the corners of the first Brillouin zone are usually organized in two sets, the Dirac points 
K and K'. They are six Dirac points in the total, but only two are worth considering due 
to the periodicity of the momenta in the Brillouin zone [2|. 

It turns out that graphene's low energy excitations are relativistic corresponding to 
massless, quasi-free fermions that can be theoretically described by the Dirac equation for 
these particles [s]. Thus, one considers the Dirac equation at the vicinity of the Dirac 
points K and K'. Expanding the dispersion relation around these points, one has to a first 
order approximation a linear relation, which, for small energies, gives origin to the so-called 
Dirac cones. These cones imply that graphene can be seen as a conventional semiconductor, 
given that there is no gap between conduction and valence bands. 

We start with Dirac equation, 

zh^ = Hn^P , (1) 

where the wave function ip in the graphene case describes the electron states around the 
Dirac points K and K', and the Dirac Hamiltonian is given by 4| 

HD = c[d-P + (3mc] , (2) 



where a and /3 are the Dirac matrices and P = {—ihdx-, —ihdy, 0). In the case of graphene, 
one has massless particles that move through the honeycomb lattice with a velocity vp ^ 
10^ ms~^, the so-called Fermi velocity. Thus, for instance, the Dirac Hamiltonian around 
the Dirac point K reads 

H = VFa- P . (3) 

where, a = {(Jx,cry,az) and ax, Uy and are the Pauli matrices. The same Hamiltonian 
can be written to the Dirac Point K' with a* given by a* = {ax, —ay,az)- We can now 
write the Hamiltonian for the two Dirac points as j4|, 

\ 



H 



D 




vf 












Px + iPy 



\ 



(4) 



Px- iPy 

Px + iPy 



Px-iPy 

As mentioned, the wave function ip consists of two components, one describing K and 
other K' . Moreover, for these two Dirac points one has an eigenvector that describes the 
probablility of an electron state to be on sub-lattice A in the upper component, or on the 
sub-lattice B in the lower component of the eigenstate. Thus, 

where ip^ and ip^' are two dimensional eigenstantes. 



(5) 



K 



1 4,-'\ 

















(6) 



To obtain the dispersion relation for the energy one has to solve the following eigenvalue 
problem for each Dirac point. 



K 



(7) 



Ek^ , 

We have for the eigenvalues, for the two Dirac points K and K', [4| 

Ek,k' = ±hvF\k\ (8) 

and to evaluate the eigenvectors of the system we consider Eq. ([7]) in the momentum 
representation, 

{a-k)tlj^{k) = ±^^{k) (9) 

{a* ■ k)^""' (k) = ±^^'(^) , (10) 
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where a and a* are defined above. Moreover, these eigenvectors in the momentum space, 
for the two Dirac points, have the form \^ 



1 



a/2 I -i-gi'/'fc/a 



K' 



1 



a/2 1 -tg-i</'fc/2 



where (pk is the polar angle of vector k in the momentum space. The signs ± correspond 
to the eigenvalues of the energy spectrum for each Dirac point, given by Eq. ([8]). One 
clearly sees that the pseudo spinor ip{k) has a definite pseudo-helicity. Furthermore, the 
existence of chiral symmetry, in theprevious eigenvectors, allows for the observation of the 
anomalous quantum Hall effect 3|-l6| . This is a striking feature of graphene, which instead 
of the usual Landau levels for a semiconductor in a magnetic field, it exhibits a different 
degeneracy factor. This is because the fundamental energy level of the graphene has two 
valleys flfl- 

It has been argued in various instances that noncommutativity should be considered 



as a fundamental feature of space-time at the Planck scale P]. Various noncommutative 



field theory models js-llll have been discussed as well as many extensions of quantum 
mechanics 12l-l21|. Of particular interest is the so-called phase-s pac e noncommutativity 
which has been investigated in the context of quantum cosmology 221], black holes physics 
and the singularity problem 23l-l25ll. The phase-space noncommutative algebra is given by 

Q 

[xi, Xj] = i9ij, [xi,Pj] = ihSij, [pi,Pj] = irjij, hj = 1, d (12) 

where rjij and 6ij are antisymmetric real constant {d x d) matrices and 6ij is the identity 
matrix. The key property of this extended algebra is that it is related to the standard 
Heisenberg-Weyl algebra: 



a; j , Xj 



0, 



ih5ij, 



P'i,p'j 



0, 



1, ...,d , 



(13) 



by a class of linear (non-canonical) transformations: 



(14) 



More recently, a new representation of noncommutativity has been proposed 261]. This 
representation uses the Pauli matrices as the fundamental elements of the noncommutative 
algebra. It allows for a noncommutative extension of graphene by treating it as a quantum 
charged particle subjected to an electromagnetic field [27]. 



Since graphene is a two-dimensional system it is expected to supply an interesting model 
where to test noncommutativity. First, because it is a special system where quantum rel- 
ativistic phenomena, typical of high-energy physics, arise at low-energies. Second because 
in non-relativistic quantum mechanics time is always a commutative variable and so non- 
commutativity, in this context, can only be applied to the spatial variables. Since any 
sympletic form in an odd dimension is always degenerate, after a linear transformation 
we can only have noncommutativity in two dimensions. Graphene provides a real two- 
dimensional physical system. Various bounds and inequalities, which appear in the context 



19|. 



of two-dimensional noncommutative quantum mechanics, where derived in 

In this work, we extend the model of graphene to a noncommutative phase-space setting. 
We consider graphene in an external magnetic field and a noncommutative geometry. In 
section II A we review the concepts behind the graphene model in the presence of a magnetic 
field. Then, in section II B we present the noncommutative extension of the graphene in 
a constant magnetic field. In section II C we examine the effect of noncommutativity 
on graphene's anomalous quantum Hall effect. In section III we compare the experimental 
values obtained for the energy levels of graphene with our theoretical predictions and obtain 
a bound for the noncommutative parameter 77. Finally, in section IV, we summarize the 
main conclusions and results. 



II. GRAPHENE IN AN EXTERNAL MAGNETIC FIELD 

In this section, we start by obtaining the energy dispersion relation for a layer of graphene 
subjected to an external constant magnetic field. We then show how to extend the problem 
to a phase-space noncommutative setting, and consider the implications of this extension 
on the anomalous quantum Hall effect. 



A. The Commutative Case 

Let us consider a layer of graphene in a external constant magnetic field, B = Be^. In 
our units c = 1. We introduce the 5-field through the minimal coupling to the vector 
potential, such that B = V x A 

P-^P-eA, (15) 
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where 



- B 



(16) 



and e is the charge of the electron. Thus, for the two Dirac points K and K' the Hamilto- 
nians read {2!: 



Hk = vf 



' Vx-Wy^^{y^^^) 

y^Px + iPy + - ix) 

/ n „ I „•„ I e_B/ 



(17) 



Hk' = Vf 



(18) 



Px + iPy + ^{y - ix) 

Px - iPy + ^{y + ix) 
and the energy eigenvalue equation for the wave function Eq. ([5]), at the Dirac point K, is 
then: 

hVF 



X 



hv} 

T 



where ^ 



\/21b is an auxiliary variable and Ib is the so-called magnetic length 



(19) 

a. 



Redefining the variables as, 



X 



(20) 



the system, Eq. (fT9l) . turns into 



hVF 



i {di -x) + {-dy + y)] 
t{d, + x) + idy + y)]<l)^ = EK<p'' . 



(a) 
ib) 



Solving Eq. (pTb) for 0^ and substituting into Eq. (l2Tb). yields 



-I (d, -x) + i-dy + y)] [-1 (d^ + x) + (dy + y)] 0^ = , 



(21) 



(22) 



which is a second order equation analogous to the energy eigenvalue equation for the quan- 
tum harmonic oscillator in two dimensions. It can be solved using the set of annihilation 
and creation operators. 



ax = -^{x + di) , 4 = -^{x - ds) , 



(23) 



However, it is more convenient to use left/ right operators, as in the case of the quantum 
harmonic oscillator in two dimensions. These new operators are responsible not only for 
adding (subtracting) a quantum of energy, but also for adding (subtracting) a quantum of 
angular momentum, h, in the direct direction (right operators with subscript d) or in the 
inverse direction (left operators with subscript e). Furthermore they simplify considerably 
the calculations. Thus, one introduces the left/ right operators as. 

So, Eq. fl22|) can now be rewritten in terms of these left / right operators, 

4(^)'(aK)0^ = i?^0^. (25) 

One sees that the Hamiltonian depends only on the left operators. Note that for the 
angular momentum one has Iz = fi{nd — rig) = km, where rig and rid are the eigenvalues 
associated to the left and right number operators {Ne = a\ae^ Nd = aj^a^) respectively; m 
is the eigenvalue associated to the angular momentum operator in the z-direction. Let us 
suppose that 0^ = 0^^^ then the energy spectrum for the Dirac point K is given by 

= ±2— ^v^=±^V^v^, (26) 

with Tie = 0, 1, 2, .... 

Finally, if we consider the eigenvalue problem, Eq. (pTb). and substitute the egeinvalue 
Ek obtained in Eq. fl2^ . we get the eigenstates for the Dirac K point. 

/ 



We also included the eigenstates for the Dirac point K' which can be obtained using the 
same method and display a dispersion relation identical to the one given by Eq. ( 126|) . 

B. The Noncommutative Case 

In this section we consider a phase-space noncommutative algebra as in Eq. f fT2|) . In 
order to relate the noncommutative variables {xi,pi) with the commutative ones {x[,p'j) we 



("e,nd) 



(ne-l,na) 
kB 



(27) 
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use the following SW map 



m 



[see also 



18|) 



(2^ 



where 6 and rj are real constant parameters. The noncommutative variables then satisfy 
the algebra 



ie. 



07] 



(29) 



2. This is the noncommutative algebra Eq.(lT2] with an effective Planck constant 
171 . Il8| . It reduces to Eq.(IT2] exactly when ^ = 0, this being the case we are going to 
consider. 

Using the same potential vector as in section IIA, Eq. f fT6|) . and substituting the non- 
commutative variables by the commutative ones, through the SW map, Eq. (!28|) . we get, 

/ 



eB , 



\Py 2 ^ 



where 



A 



1 



V 



Ah ) ' V eBh 
Thus, for the Dirac point one gets the following Hamiltonian: 

/ 



(30) 



Hk = vf 







Kp'. + ¥y) + ff^iy' - 



Kp'.-¥y) + fKy' + 



IX 







(31) 



A straightforward comparison with Hamiltonian Eq. (fT7|) shows that noncommutativity 
reveals itself through constants A and /i. For future convenience one introduces the constant 
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I2h\_^ 2\ 



(32) 



In what follows we shall consider only momenta noncommutativity, since in the general 
Dirac problem, configuration space noncommutativity leads to the breaking of gauge sym- 
metry [ll|, a symmetry preserved in the graphene lattice [29] . Indeed, if one evaluates the 
velocity of a charged particle. 



V = -[Hnc^A 



(33) 



we obtain an extra term depending on the 9 parameter, and not the expected result v 
vf(T- Hence A = 1 and 



2h 1 , 2 



(34) 



Thus, the equations to be solved are now. 



7 
7 



-i 7<9x' 



X 



1, 



X 



-idy 



7. 



-i 7(9^. + - 
V 7 



E 



(35) 



Following the strategy discussed in the last section, one obtains for the energy spectrum 



E 



NC 
K 



V 



± 



Ib \ 



(36) 



where rif. is a non-negative integer. We clearly see that the energy at the Dirac point 
K depends explicitly on the noncommutative parameter associated with the momenta. 
Moreover, one concludes that the noncommutative effect is coupled with the magnetic 
field. Of course, the same dispersion relation for energy is found to the other Dirac point, 
K' . Furthermore, the eigenvectors for the two Dirac points can be evaluated. One has 



/ 



lA 
B 

ne-l,na) 



(37) 



The eigenvectors can be written in polar coordinates 28 1 



{31 



where L = A, B , a > is some constant, = + y^, J^{p) is some polynomial of p and 
m is the quantum number associated to the angular momentum l^, which is the integer 
{m = nd- He). 

Clearly, the zero-energy level for this system is E^'-' = E^f = 0. However, the eigen- 
vectors for the fundamental level are 

/ 



K 



no) 



( A.A 











(39) 



meaning that the zero energy for the two Dirac points are associated with two linearly 
independent electron states. 
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C. Quantum Hall effect and noncommutativity 



We conclude from the discussion in the previous section that an external magnetic field 
perpendicular to a graphene sheet renders the energy spectrum discrete. It is known that 
a charged particle subjected to an electromagnetic field has a discrete energy spectrum, 
the Landau quantization. However, graphene is a special case, since it is actually a Dirac 
problem instead of a Schrodinger one ^. The difference between each Landau level is 
considerable. This difference in the energy levels has implications for the quantum Hall 
effect (QHE). 

The QHE is observable in two-dimensional metals, as for instance in bound low- 
temperature surfaces where electrons are constrained to two dimensions That occurs 
when the temperature is drastically reduced, the Hall resistivity becomes independent of 
the magnetic field, and a quantized Hall plateau is formed. Experimentally, the Hall con- 
ductivity is given by 



(Txy = -j^i^ , (40) 
where u is an integer number and one has the Integer quantum Hall effect (IQHE). However, 
the QHE observed in graphene is anomalous (AQHE) js, O]. The difference between the 
IQHE and the AQHE lies in the fundamental energy level, i.e for the IQHE the first Landau 
level is observable at zero energy, but for AQHE the first Hall plateau appears when the 
lowest Landau level is half filled, and the conductivity takes the form, 

axy = ±A[iy+-jj , (41) 

where the factor 4(z/ +1/2) is evaluated by taking in account the presence of a zero mode 
shared by two Dirac points, an that there are 4(z/ + l/2) occupied states that are transferred 
from one edge to another. Notice that this effect shows up at room temperature, contrasting 

nn 

with the usual QHE in semiconductors which is typically a low temperature effect, [2|, |3[. 
Following Ref. ^|, whenever the Fermi level lies in a gap, the Hall conductivity is given 

by 

— (42) 

where n^ep) is the density of states, and ep is the energy levels of graphene given by Eq. 
dH]). Thus, in the usual commutative case, the energy levels are given by Eq. ( 126|1 and so 
the density of states is 



eB , ^ eB 

ep = ± nvp\l2—ne =^ n[ep) = — . (43) 
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The same strategy can be used to evaluate the Hall conductivity for the noncommutative 
case. In this case, the energy levels are given by Eq. (36), where the noncommutative 
correction is explicit. In what concerns the density of states, it seems rather logical that it 
is affected by the noncommutativity and the simplest way to incorporate this dependence 
is through the expression: 

ni^F) = ^ fl + . (44) 



h V eBh 
Thus, 

/e / T] \ eB r]eh \ , 

where (7^ = 4(z/ + 1/2) is the degeneracy factor. One concludes that even though the 
momenta noncommutativity induces a change in the energy levels of the graphene electrons, 
it does not affect the Hall conductivity. 



III. COMPARISON WITH EXPERIMENTAL RESULTS 



In this section, a bound on the noncommutative parameter, rj, can be obtained using 



the available experimental results. In Ref. 



3l|, infrared (IR) spectroscopy in the presence 



of a magnetic field was used to resolve the levels of the Landau spectrum for one single 
layer graphene. Two resonances were resolved for magnetic fields up to -B = 18 T ^, 
and their energy position was shown to scale as \/B with a slope corresponding to a 
vp = (1.12±0.02) X 10^ ms~^ for a particular energy Isij. This value for the Fermi velocity 
is related with the transition from n = — 1 to n = (in the case of holes) and n = to 
n = 1 (in the case of electrons). These transitions correspond to the filling factor u = —2 
in the IQHE, the lowest Landau level transitions possible in graphene. 
The energy for the Landau levels n = 1 or n = —1, 

-hv 



Ek = ±V2—^ = ±J2ehvj,B , (46) 
where + is for electrons and n = 1, and — is for holes and n = —1. For these levels one has 



Ek = ±(172 ± 3) meV 3l|- Thus, if the noncommutative energy spectrum for graphene 



is given by Eq. f l36|) . and considering that the uncertainty in the energy is at most 6 meV, 



^ Notice that for magnetic field above 10 T, the Zeeman energy g^sB is negligible in comparison with the 
energy of the Landau levels Q] . 
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one can obtain a bound for the noncommutative parameter rj. Using Eq. f l36p it follows 
that: 

r7-| < 0.069 . (47) 
Thus, the noncommutative parameter t] satisfies 

r/ < 2.1 X 10"^^ kgVs-"^ 

^ 07 < 8.6 eV/c . (48) 
Naturally, this bound is not as stringent as the one arising from the hyperfine transition 



in the hydrogen atom, ^ < 2.26 ^eV/c, HI], one of the most accurate experimental results 
in the whole of physics. Despite of that, the above reasonings show that noncommutative 
effects are consistent with what is known about graphene physics. 

In what concerns other bounds for the momentum noncommutative parameter, notice 
that the one arising .0. *e ,.av.a«o„a> ,ua.u,n well □ a.d .o. equivalence 



principle 32| depend on an assumption about the configuration space noncommutative 



parameter, 9, and cannot the compared with the above bound without fixing a value for 6. 

IV. CONCLUSIONS 

In this work the phase-space noncommutative extension of the graphene in the pres- 
ence of an external constant magnetic field was examined. More precisely, only momenta 
noncommutativity was considered since the noncommutativity associated with the config- 
uration variables implies the breaking of gauge invariance. The introduction of momenta 
noncommutativity determines a correction of the energy spectrum of graphene in the pres- 
ence of a magnetic field. Moreover, it was shown that this noncommutativity does not 
affect the graphene anomalous quantum Hall effect. 

Finally, comparison with experimental data reveals that ^ < 8.6 eV/c, a bound that is 
not very stringent, but that indicates that there is no contradiction between noncommuta- 
tive effects and graphene's physics. These results show that momentum noncommutativity 
yields interesting results also at low-energies and that its implications are not restricted to 



quantum comology |22l] and black holes physics [23 
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